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O . Abstract. Physics potential of future measurements of atmospheric neutrinos is explored. 

' Observation of Arn^i driven sub-dominant effects and #13 driven large matter effects in 

^0 , atmospheric neutrinos can be used to study the deviation of #23 from maximality and its octant. 

Neutrino mass hierarchy can be determined extremely well due to the large matter effects. New 
physics can be constrained both in standard atmospheric neutrino experiments as well as in 
future neutrino telescopes. 
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1. Introduction 

Atmospheric neutrinos observed in the Super-Kamiokande (SK) experiment provided the first 
\& • unambiguous signal for neutrino flavor oscillations lj. The observed zenith angle and energy 

dependent depletion of atmospheric v^/v^ in SK can be explained only by v il -v T oscillations with 
Am^ = 2.1 x 10 -3 eV 2 and almost maximal mixing angle, sin 2 IO23 = 1- The results of the SK 
experiment was subsequently corroborated by the MACRO and Soudan-2 atmospheric neutrino 
Q_i| experiments and more recently by the K2K and MINOS long baseline (LBL) experiments. In 

this talk we will expound the physics potential of future atmospheric neutrino experiments using 
larger and better detectors. 
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^ ■ 2. Confirming oscillations of atmospheric neutrinos 

The "smoking gun" signal for v il -v T oscillations is the observation of the characteristic "dip" in 
L/E, predicted by neutrino flavor mixing. Although the analysis of the L/E binned atmospheric 
neutrino data in SK has been found to support the oscillation hypothesis 0, it would be 
worthwhile to make an unambiguous check using a detector with better E and L resolution. This 
can be done in large magnetized detectors, such as the proposed ICAL detector at the India- 
based Neutrino Observatory (INO) Analysis of results obtained from detailed simulations 
by the INO collaboration show that oscillations can be confirmed with a significant C.L. with 
just 250 kTy data. 



3. Precision measurement of Am^ and sin 2 623 

Both Am 2 ! and sin 2 2#23 are expected to be measured very accurately by the forthcoming LBL 
experiments. A statistical analysis of the combined data set with five years of running of MINOS, 
ICARUS, OPERA, T2K and NO^A each, reveals that Am^ and sin 2 #23 could be measured 
with a spread of 4.5% and 20% respectively at 3a [I]- The future prospective data from water 
Cerenkov atmospheric neutrino experiments with a statistics 20 times the current SK statistics 
could measure Am^ and sin 2 623 with a spread of ~ 17% and ~ 24% respectively fi? a . A large 



magnetized iron calorimetric detector such as the proposed INO detector ICAL [2], could use 
atmospheric neutrinos to measure Am^ and sin 2 823 within 10% and 30% respectively at 3a 
with a statistics of 250 kTy pj. 



4. Atmospheric neutrino experiments: Subdominant effects 

The effect of the sub-dominant terms in the Super-Kamiokande (SK) atmospheric neutrino data 
is not yet at the statistically significant level. However, the sub-dominant terms, if observed in 
a future high statistics atmospheric neutrino experiment, can be used to give information on: 

• Deviation of #23 from its maximal value 

• Octant of #23 

• sgn(Am^ 1 ) 

Assuming a constant density for the earth matter, the excess of electron type events in a water 
Cerenkov experiment such as SK is given by [SI [7] 

^ - 1 - sm 2 20% sin 2 { ^f^ j x (r cos 2 823 - 1) (1) 
+ sin 2 28fi sin 2 ( ^^ ) x (r sin 2 23 - 1) (2) 



+ sin 623 cos 023 r Re 
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* 3 A 12 exp(-i<5op) , (3) 



where L is the baseline, E is the energy of the neutrino, r = N e /N^, N e and being the 
number of e and [i events respectively in the detector in absence of oscillations and 9^, #13, 
(Amj])"' and (Am^)^ are the mixing angle and mass squared differences in matter. 

(i) The first term in Eq. @ is the Am^ driven oscillation term - which is obviously more 
important for the sub-GeV neutrino sample. Since r ~ 0.5 in the sub-GeV regime, this 
term brings an excess (depletion) of sub-GeV electron events if 823 < ir/4 (623 > tt/4). It 
can thus be used to study the maximality and octant of 823 through the sub-GeV electron 
sample [H1IE]- 

(ii) The second term is the #13 driven oscillation term. Being dependent on sin 2 823 , this term 
goes in the opposite direction to the first term. Therefore for sub-GeV neutrinos, larger 
#13 would imply that the effect of the first term would get negated by this term. However 
for multi-GeV neutrinos, there will be large matter effects inside the earth and this term 
dictates the electron excess. The sin 2 823 dependence of this term could then be used to 
study the maximality and octant of 823 through the multi-GeV electron sample (see also 
[H]). Since matter effects bring in sensitivity to the sgn^Am^), this term can be used to 
study the mass hierarchy. 

(iii) The last term is the "interference" term [7j, which depends on 5c p- The effect of this term 
could be to dilute the effect of the first two terms and spoil the sensitivity of the experiment. 
However, being directly dependent on 5cp, this term also brings in some sensitivity to the 
CP phase itself (HHg. 

The depletion of the muon events in the limit of Am^ = is given by 1 

1 ~W = (Pl» + Pl») + (Pl»)'^ 823(^823--) , (4) 

1 The approximation of taking a vanishing Am^ has been made in Eq. 10 only for the sake of simplicity, since 
the main subdominant effect in the muon neutrino channel comes from earth matter effects, which are large for 
multi-GeV neutrinos for which AmJi dependence is less importance. The results presented in the later sections 
have been obtained using the full numerical solution of the three-generation equation of the atmospheric neutrinos. 



P},„ = sin 2 6fl sin 2 26» 23 sin 2 
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(P^Y = sin 2 20% sin 2 ^f 1 , (7) 

where ^4 = 2\^2GpN e E is the matter potential. For very small values of #13, there is very little 
matter effect and we can see that P 2 ^ is the dominant term in the survival probability. Since 
this term depends on sin 2 2#23 we do not expect octant sensitivity in absence of matter effects 
from experiments probing the Pun channel alone. However, if #13 is not too small, neutrinos 
which travel through large baselines and hence large matter densities inside the earth, undergo 
large matter effects. The mixing angle 6^ increases in matter and the third term {Pan)' becomes 
important as well. Since this term has a strong dependence on sin 2 #23 rather than sin 2$23 j we 
expect the Pnn channel to develop sensitivity to the octant of #23 in presence of large matter 
effects p]] • Probing matter effects in the resultant muon signal in the detector will also provide 
us with information on the neutrino mass hierarchy 11 2 j . 

Since matter effects are large for higher energy neutrinos, we expect that multi-GeV 
atmospheric v^jv^ events can be used for this purpose. However, unlike in the case for matter 
effects in the P^ e channel, both the magnitude and sign of the earth matter effects in the Pnn 
channel depends crucially on L and E. The largest effect of earth matter comes for neutrino 
travelling L ~ 7000 km with E ~ 5 GeV. The matter effects changes sign rapidly with L and E 
- with A(P„„) < and A(P W ) > at the maximum and minimum respectively of Pnn- Thus, 
in order to see the matter effects one needs to bin the data judiciously both in energy and zenith 
angle. 

Very good energy and zenith angle detector resolution is expected for the magnetized iron 
calorimeters. Therefore, fine binning would allow such detectors to observe matter effects in the 
muon signal. Since large matter effects appear only in either the neutrino or the antineutrino 
channel, the magnetic field which allows for charge discrimination, further helps these type of 
detectors to observe earth matter effects in the muon channel. However, unless the iron plates 
of the detector are thin enough, it would not be possible to detect electrons in these type of 
detectors. The current INO-ICAL design does not allow for it and therefore would observe muon 
events only. Another restriction for these detectors come from the relatively higher threshold, 
which allows for the detection of only multi-GeV v^jv^. 

Water Cerenkov detectors have the advantage that sub-GeV neutrinos can be detected. 
However, the energy resolution is worse than that for iron calorimeters. For the results shown 
here, the data is binned in sub-GeV and multi-GeV bins and therefore the matter effects in 
the Pnn channel get largely averaged out. This means that one would see very small residual 
matter effects in the multi-GeV muon sample. However, matter effects in the P^ e channel do 
not change sign over most of the relevant range of E and L in the multi-GeV regime. Therefore, 
multi-GeV electron sample has large matter effects and can be used to study the deviation of 
#23 from maximality and its octant as well as the neutrino mass hierarchy. 



5. Is the mixing angle #23 maximal? 

The measurement of both the magnitude and sign of the deviation of sin 2 #23 from its maximal 
value 0.5 is of utmost theoretical importance. To quantify the deviation of the true value 
of #23 from its maximal value, we introduce the function D = | — sin 2 #23- The magnitude 
\D\ gives the deviation of sin 2 #23 from its maximal value, while sgn(D) gives the octant of 
sin 2 #23- The best current limit on \D\ comes from the SK atmospheric neutrino experiment 
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Figure 1. The regions of Amg^true) and sin 2 #23 (true) where maximal #23 mixing can be 
rejected at la (inner bands), 2a (middle bands) and 3d (outer bands) C.L. The left-hand panel 
shows the sensitivity expected from the combined data from the LBL experiments. The 
middle panel 5 a shows the sensitivity expected with atmospheric neutrinos in a megaton water 
detector (SK50). The extreme right-hand panel |T(]] shows the corresponding reach expected 
from 500 kTy atmospheric neutrino data in large magnetized iron detectors. The true value of 
$13 is assumed to be zero. 
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Figure 2. Left-hand panel same as in Fig. ^ Middle panel (for water detector) and right-hand 
panel (for magnetized iron detector) have been drawn assuming that sin 2 #13 (true) = 0.04. 



giving \D\ < 0.16 at 3a while sgn{D) is almost unknown at present. Fig. ^ shows the 
potential of atmospheric neutrino experiments to test the deviation of #23 from maximality and 
compares it with the reach of the combined data from the current and next generation long 
baseline experiments. The combined long baseline data include five years of running each of the 
MINOS, ICARUS, OPERA, T2K and NOz/A experiments. The middle panel gives the sensitivity 
to \D\ of atmospheric neutrino experiments with water detectors for a 4.6 Megaton-yr statistics, 
while the left panel shows the corresponding sensitivity of atmospheric neutrino data in large 
magnetized iron detectors with 500 kTy statistics. At Am 2 1 (true)= 2.5 x 10 -3 eV 2 , it should 
be possible to measure \D\ within 19% and 25% at 3a with atmospheric neutrinos in water and 
iron detectors respectively. This is slightly weaker than the sensitivity of the combined long 
baseline experiments, where it should be possible to measure \D\ within 14% at 3a. However, 
note that all the results presented in Fig. ^have been obtained assuming that the true value 
of #13 was zero. Results for atmospheric neutrino experiments when the true value of #13 is not 
zero is shown in Fig. [21 For non-zero #13, presence of earth matter effects in the Pau channel 
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Figure 3. Plot showing the octant sensitivity as a function of sin 2 #23 (true), for an atmospheric 
neutrino experiment with large magnetized iron calorimeter (left-hand panel) and megaton water 
detector (right-hand panel). 
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Table 1. A comparison of the potential of different experiments to rule out the wrong #23 
octant at 3a (1 dof). The third column gives the condition on the true value of sin 2 #13 needed 
for the 623 octant resolution. 



brings in a marginal improvement in the sensitivity of atmospheric neutrino experiment with the 
magnetized iron detector. For the megaton water atmospheric neutrino experiment, very large 
earth matter effects in the P^ e channel bring in significant improvement in the determination of 
|D|, making this experiment comparable/better than the long baseline experiments for studying 
the deviation of #23 from maximality. 

6. Resolving the #23 Octant Ambiguity 

If the true value of #23 is not 45°, then the question arises whether #23 > (D positive) or < 7r/4 
(D negative). This leads to an additional two- fold degeneracy in the measurement of the mixing 
angle #13 and the CP phase 5c p in LBL experiments. This ambiguity is generally regarded as 
the most difficult to resolve. As discussed before, the presence of earth matter effects in the 
zenith angle and energy binned atmospheric v^jv^ data in magnetized iron detectors opens up 
the possibility of probing the octant of #23 [TO] - On the other hand atmospheric v e /i> e data 
in water detectors could also give information on the octant of 023, both through the Am^ 
dependent subdominant term in the sub-GeV sample [f3 Ej, as well as through earth matter 
effect in the multi-GeV events, as discussed above. One could hence combine the atmospheric 
neutrino data, in either megaton water detectors [HI E] > or in large magnetized iron calorimeters 
with data from long baseline experiments to resolve parameter degeneracies. 

In order to obtain the limiting value of sin 2 #23 (true) which could still allow for the 
determination of sgn(D) we define 



Ax 2 = X 2 (sin 2 023 (true) , sin 2 #13 (true) , others(true) ) — X 2 (sin 2 6*23 (false) , sin 2 #13 , others) , (8) 



with sin 2 #23 (false) restricted to the wrong octant and 'others' comprising Am^, Am|i, sin 2 #12 
and 5c p- These, along with sin 2 #13 as well as sin 2 #23 (false), are allowed to vary freely in 
the fit. Fig. El shows the results of a statistical analysis based on simulated data from 
atmospheric neutrinos with 500 kiloton-yr exposure in a large magnetized iron calorimeter (left- 
hand panel) and 4.6 Megaton-yr exposure in a water Cerenkov experiment (right-hand panel). 
For large magnetized iron detector we show results for four different values of sin 2 #13 (true) , 
assuming a normal mass ordering. For a given sin 2 #13 (true), the range of sin 2 #23 (true), for 
which sin 2 023 (false) can be ruled out with atmospheric neutrinos in magnetized iron detector 
is given in Table ^ This can be compared to the sensitivity possible with water Cerenkov 
detectors, shown for a true normal hierarchy in the right-hand panel of Fig. |3] and Table ^ 
where octant determination can be done reasonably well even if sin 2 #13 (true) was zero (Sj. 
However, if sin 2 #13 (true) was non- vanishing and reasonably large, the octant sensitivity of this 
experiments gets significantly boosted through earth matter effects appearing in the multi-GeV 
electron sample. 

7. Resolving the Ambiguity in Neutrino Mass Hierarchy 

Large earth matter effects in atmospheric neutrinos can be exploited to probe the sign of Am^. 
Fig. 0] 12 shows the sensitivity to sgn^Am^) expected in a magnetized iron calorimeter, 
with 4000 observed upward going events. The data corresponds to a normal (solid lines) and 
inverted (dashed lines) hierarchy and the curves show the x 2 and hence the C.L. with which the 
wrong hierarchy can be ruled out. The red lines correspond to an analysis method where all 
the oscillation parameters are fixed in the fit. The blue lines show the results of the fit where 
external priors for the oscillation parameters have been used. The green lines correspond to the 
case where all oscillation parameters are allowed to vary freely in the fit. The left-hand panel 
is for muon events in a detector with 15% energy and 15° zenith angle resolution, the middle 
panel is for muon events with 5% energy and 5° zenith angle resolution, while the right-hand 
panel is for electron events. For vanishing #13 the matter effects vanish giving x 2 = 0- As 
#13 increases, matter effects increase, thereby increasing the sensitivity of the experiment to 
hierarchy determination. For a INO-ICAL like detector, where energy resolution is expected 
to be around 15% and zenith angle resolution of about 15°, the wrong hierarchy can be ruled 
out at 2a using the muon events, if sin 2 2#i3(true)= 0.1 and sin 2 #23(t r ue)= 0.5, and where the 
information from the other long baseline experiments on the oscillation parameters have been 
included through the priors. Comparison of the left-hand with the middle panel shows that the 
sensitivity to hierarchy increases if the detector resolution is improved. Comparison of the left- 
hand with the right-hand panel shows that the sensitivity to hierarchy increases if the detector 
could detect electron type events as well. And of course since matter effects increase with 623, 
the sensitivity to hierarchy increases as the true value of 623 increases. 

The sgn(Am 2 1 ) can be done using the excess in the multi-GeV electron sample due to 
earth matter effects in water Cerenkov detectors. The wrong hierarchy can be ruled by a 4.6 
Megaton-yr data in such an experiment at more than the 2a limit if sin 2 2^13 (true) = 0.1 and 
sin 2 #23 (true) = 0.5 (see also (Sj)- This is comparable to the sensitivity of the magnetized iron 
detectors as discussed above. However, since water detectors use the excess in electron events for 
multi-GeV neutrinos, which in turn have large matter effects in the P Me channel, they therefore 
depend also on the CP phase 5cp as discussed before. If the value of 5cp is allowed to vary 
freely in the fit then the sensitivity gets affected and decreases appreciably. 

8. Looking for new physics 

There are a variety of new physics scenarios which could manifest themselves as subdominant 
effects in oscillations of atmospheric neutrinos. Each one of these have a distinctive L/E 
behavior, while oscillations go linearly with L/E. This can be probed directly using the L/E 
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Figure 4. Ay 2 for the wrong hierarchy as a function of sin 2 2#i3(true). See |12j for details. 



binned data in large iron detectors ° r by comparing the low energy contained events with 
the high energy upward going muons in large water detectors • Atmospheric neutrino events 
which constitute a background for the neutrino telescopes such as IceCube, can also be used every 
effectively to constrain new physics. Neutrino telescopes look for neutrinos in the (10 _1 — 10 ) 
TeV range, for which standard oscillations of atmospheric neutrinos are negligible. Hence, any 
E or L dependence in the data would signify new physics. Likewise, the absence of any E and 
L dependence can be used to constrain the new physics parameters |17j . 

9. Conclusions 

In this talk we explored the physics potential of future measurements of atmospheric neutrinos. 
Observation of Am^i driven sub-dominant effects and #13 driven large matter effects in 
atmospheric neutrinos can be used to study the deviation of #23 from maximality and its octant. 
Neutrino mass hierarchy can be determined extremely well due to the large matter effects. New 
physics can be constrained both in standard atmospheric neutrino experiments as well as in 
future neutrino telescopes. 
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